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(57) Abstract: A system for monitoring mechanical waves from a machine which in operation has moving particulate matter therein, 
the system including at least one sensor located on the machine at a location away from the central axis of the machine, the sensors 
being for sensing acoustic waves and including a transmitter for transmitting signals representing the sensed mechanical waves to 
a receiver at a location remote from the sensor(s), a data processor connected to the receiver for receiving signals from the receiver 
which signals represent the mechanical waves and processing the signals to produce output signals for display on a display means, 
wherein the output signals for display represent one or more parameters indicative of mechanical waves emitted from the machine 
over a predetermined period of time. 
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A SYSTEM FOR MONITORING MECHANICAL WAVES FROM A MOVING 
MACHINE 

Field of the Invention 

This invention relates to the monitoring and 
• 5 analysis of surface vibration waves generated by the 
operation of material processing equipment. 'The invention 
is of particular application to the non- intrusive 
monitoring and control of mineral processing equipment, 
such as tumbling mills, where the equipment component being 

10 monitored is in motion. The invention also has application 
to equipment with stationary components but moving mineral 
particles or pulp flows, such as crushers and 
hydrocyclones . The invention was initially developed for 
condition monitoring and process control of Semi -Autogenous 

15 Grinding (SAG) mills. 

Background of the Invention 

Acoustic emissions and surface vibrations 
monitoring, have previously been used to investigate and 
control t^he performance of mineral processing unit 

20 operations. Control of power draft in autogenous grinding 
(AG) mills and SAG mills has traditionally been via load 
cells estimating the charge mass. However, acoustic 
emissions from dual microphone systems have been used to 
monitor the changing level of impact of the charge on an AG 

2 5 mill shell. A pair of microphones were mounted at 

approximately 30° from the bottom and 30° from the centre 
line of the mill. The position of the microphones was such 
that the upper unit was above the normal level of charge 
impact on the liner while the lower unit was below the 

3 0 same. If the load level rises in the mill, the point of 

impact moves toward the upper microphone and away from the 
lower. If the load level drops, the converse applies. 
Therefore changes in load level are registered by variation 
in the comparative sound intensity at the two microphones. 
3 5 The resulting estimate of the load volume is correlated 
with the power draft and used to control the feed rate in 
order., .to maintain optimal milling conditions and maximum 
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power draft. It was shown that microphone signals are much 
more sensitive to load change than the load cell. However, 
the method is crude in that it uses sound intensity at only 
two fixed points outside the mill. The intensity of sound 
5 at these two positions may be considerably influenced by 
events both outside the mill and at a variety of locations 
within the mill. The technique therefore only permits 
qualitative investigation of the state of the charge inside 
the mill. 

10 Acoustic emissions are also known to be 

indicators of pulp density and viscosity. The dual 
microphone study of AG mill acoustic emissions showed that 
the sound intensity emanating from the charge region (lower 
microphone) was correlated with the pulp density- The lower 

15 microphone sound intensity was used to control water 
addition rate. Low pulp density was thought to result in 
higher transmission of noise and increased media/media and 
media/liner collision events. Meanwhile at higher pulp 
density grinding action was thought to be inhibited by the 

2 0 increased pulp viscosity, reflected in lower noise 

intensity. Estimation of effective pulp density and 
viscosity via the magnitude of acoustic emissions has also 
been achieved for laboratory batch ball mills. Results 
suggest that changes in mill noise can be used to identify 
25 the pulp rheological regime and potentially used to 
optimise grinding efficiency. Mill sound noise has also 
been shown to indicate charge size distribution, ore 
breakage rates, and ore character in batch ball mills 
(Watson, 1985; Watson and Morrison, 1985). 

3 0 Acoustic emissions monitoring has also been used 

to analyse hydrocyclone performance. A shear structure 
piezoelectric type acoustic sensor was mounted halfway 
along the conical section of a 5" hydrocyclone body. The 
digitised signal was sampled at 200 0 Hz and a Fast Fourier 
35 Transform (FFT) algorithm used to derive the Power Spectral 
Density function (PSD) for analysis of acoustic emission 
characteristics in the frequency domain. Features of 
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- 3 - 

acoustic emissions were analysed for varying feed solids 
concentrations and pressure. Results indicated significant 
spectral features in the frequency range from DC to about 
50 Hz and between 3 0 and 45 Hz. The height of these 

1 5 spectral features was sensitive to operating conditions. It 
was conjectured that the spectral structure is related to 
features of the hydrodynamics inside the hydrocyclone local 
to the sensor. A stepwise regression analysis technique was 
used to derive linear relationships between the operating 

10 parameters of the cyclone and the spectral and statistical 
characteristics of the acoustic emissions. The signal 
measures used in this analysis were for the time domain 
maxima, mean, standard deviation, rootmean- square, skewness 
and kurtosis, and for the frequency domain the first 52 

15 spectral components of the PSI) . The model was then used 
for reasonable predictions of hydrocyclone feed pressure, 
solids concentration, mass and volume flow rates and 
underflow concentration. This investigation showed that 
non-invasiye acoustic emission measurement coupled with 

20 multivariate statistical analysis techniques are a useful 
tool for monitoring the bulk characteristics of both 
process and equipment, in this case hydrocyclone operation. 

Vibration monitoring and signal analysis have 
been used to study the feed distribution characteristics of 

25 parallel Dense Medium (DM) cyclones in a coal preparation 
plant. The method is based on the concept that the 
monitoring of vibrations on the external surface of the 
cyclone can yield the frequency and strength of particle 
impacts (particularly for larger particles near entry and 

3 0 exit points) . Accelerometers for measuring vibrational 
accelerations were mounted near the feed inlet, underflow 
spigot and overflow cap. Relatively large vibrations were 
noted in the region of the overflow cap, reflecting the 
energy of particle/wall impacts in that region due to the 

35 flow regime within the cyclone. Results indicated that 
vibration measurements are a superposition of a large 
number of transients caused by individual particle impacts. 
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Summary of the Invention 

According to the present invention a system is 
provided which is capable of monitoring acoustic emissions 
from a moving machine, 
5 The invention includes within its scope systems 

for monitoring moving substances within a stationery or 
moving machine. 

The system in broad terms is directed at 
monitoring mechanical wave emissions from inside the 
10 machine and the surface of the machine , as well as 
associated components of the machine which are affected by 
the machine's operation. 
Summary of the Invention 

According to a first aspect, the invention 
15 provides a system for monitoring mechanical waves from a 
machine which in operation has moving particulate matter 
therein, the system including at least one sensor located 
on the machine at a location away from the central axis of 
the machine, the sensor (s) being for sensing mechanical 

2 0 waves and including a transmitter for transmitting signals 

representing the sensed mechanical waves to a receiver at a 
location remote from the sensor (s), a data processor 
connected to the receiver for receiving signals from the 
receiver which signals represent the mechanical waves and 

25 processing the signals to produce output signals for 
display on a display means, wherein the output signals for 
display represent one or more parameters indicative of 
mechanical waves emitted from the machine over a 
predetermined period of time. 

30 It is preferred that the output signals represent 

a number of acoustic events occurring within the machine, 
amplitudes of the acoustic events and data relating to the 
position of the acoustic events. 

It is preferred that the system includes a 

3 5 plurality of sensors each for detecting acoustic emissions 

from inside the machine and from the surface of the 
machine . 
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It is preferred that the system includes a 
plurality of sensors spaced around the periphery of the 
machine to enable polar co-ordinates of the origin of 
emissions to be located. 
5 It is preferred that the sensors are eguispaced 

around the periphery of the machine. 

According to one embodiment of the invention 
sensors are arranged in an array around the machine and 
along the length of the machine to enable a three- 
10 dimensional co-ordinate axis to be plotted of the location 
of the origin of emissions from the machine. 

It is preferred that these sensors are removably 
attached to the outer surface of the machine. 

According to one embodiment the transmitter 
15 associated with each sensor is located away from a detector 
part of the sensor. 

According to one embodiment one or more 
transmitters of the sensors are removably attachable to the 
machine . 

20 It is preferred that the system includes at least 

one proximity detector for monitoring the location of the 
sensors at a predetermined time. 

According to one variation of the present 
invention the data processor includes a timing means for 

2 5 calculating the location of the sensor (s) at a 

predetermined time. 

It is preferred that the timing means output data 
relating to the position of the sensor (s) at a particular 
time, based on data received either from the proximity 

3 0 detector (s) or/and data received relating to the movement 

of the machine. 

It is preferred that the data processor receives 
signals from the receiver, which signals include data 
relating to the frequency of vibrational events occurring 
3 5 within the machine and the amplitude of the vibrational 
events at particular locations within the machine. 

It is preferred that the sensor (s) includes an 
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accelerometer . 

It is preferred that the sensor includes a power 

supply. 

The power supply may be a solar cell. 

• 5 According to another embodiment of the present 

invention the power supply is a 12 volt gel cell 
rechargeable battery. 

Battery charging may be achieved using two solar 
panels mounted on opposite sides of the mill or 

10 alternatively other recharging methods include an inertial 
generator, offtake from the electric drive, or some other 
source that provides a continuous power supply. 

According to another aspect of the present 
invention there is provided a method of analysing 

15 operational parameters of a machine having a moving 
particulate material therein, the method including the 
steps of recording data representing a number of mechanical 
events occurring within the machine over a predetermined 
period of x time, the amplitude of the mechanical events 

2 0 occurring over the predetermined period of time and 

positional data relating to the position of the mechanical 
events occurring within the machine, displaying a graphical 
representation of the recorded data, the graphical 
representation including parameters relating to the number 
25 of mechanical events, the amplitude of mechanical events 
and the position of mechanical events occurring within the 
machine during the machines operation. 

Preferably the mechanical events include 
mechanical events. 

3 0 It is preferred that the graphical representation 

includes data on the radial and angular position of each 
vibrational event. 

Preferably the graphical representation shows the 
number of vibrational events on a polar co-ordinate graph. 
35 According to one embodiment the step of 

displaying includes displaying a histogram of variables 
relating to the number of vibrational events occurring over 
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the predetermined time. 

The present invention also includes within its 
scope graphical representation of the Recorded data in a 
rectilinear co-ordinate system and three dimensional co- 
5 ordinate system. 

According to one embodiment of the present 
invention the step of displaying includes displaying a 
graphical representation of the recorded data on a data 
processing monitor using 3D graphics which simulate the 
10 machine and provide a graphical representation of the 
moving particulate material within the machine. 

It is preferred that the step of displaying 
includes highlighting regions within the machine which have 
vibrational events liable to cause maximum damage to a 
15 component of the machine. 

According to another variation of the present 
invention the graphical representation includes a colour 
scheme for colouring regions of the graphical 
representation according to features of vibrational events 

2 0 occurring at positions of the vibrational events. 

According to another aspect of the present 
invention there is provided a method of controlling 
operational parameters of a machine having a moveable 
substance therein, the method including the steps of 
25 recording data representing a number of mechanical events 
occurring within a machine over a predetermined period of 
time, amplitude and/or frequency of the mechanical events 
occurring over the predetermined period of time and 
position data relating to the position of the mechanical 

3 0 events over the predetermined period of time, determining 

zones within the machine which are subject to predetermined 
levels of wear and altering the machine operational 
characteristics to reduce the levels of wear for the zones. 

Preferably the mechanical events includes 
35 vibrational events. 

It is preferred that the step of determining 
zones, includes processing the recorded data with data 
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relating to the substance or substances within the machine, 
and dynamic properties of the substance (s) to produce a 
level of wear indication parameter for a plurality of zones 
within the machine. 

•5 It is preferred that the method includes 

operating a data processor which includes recorded model 
data relating to wear characteristics of the machine as a 
function of a plurality of parameters which may include one 
or more combinations of the following: 

10 The number of vibrational events occurring within 

the machine, the amplitude of the vibrational events 
occurring within the machine, the position of the 
vibrational events occurring within the machine, the mass 
of the particulate material and other material substance (s) 

15 within the machine, the size of the particulate matter 
within the machine, the volume of particulate matter within 
the machine, the volume of space within the machine, the 
shape of the machine and other parameters which are likely 
to affect tyhe wear characteristics of the machine. 

2 0 Alternatively the method of controlling the 

machine includes determining the efficiency of operation of 
the operational characteristics of the machine. 

It is preferred that the method includes the step 
of processing the recorded data with other data relating to 
25 characteristics of the particulate matter and machine to 
determine the efficiency of operating characteristics of 
the machine. 

The operating characteristics could include the 
efficiency of a crushing operation over a predetermined 

3 0 period of time. 

It is preferred that the method includes the step 
of increasing or decreasing the speed of operation of the 
machine, including rotation or reciprocating motion of the 
machine and/or alteration of rate of feed of particulate 
3 5 matter to the inside of the machine. 

It is preferred that the method includes the step 
of .maximising _a .predetermined operational parameter of the 
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machine. This may include maximising the amount of 
crushing of a material within a machine that is a crusher. 
Preferably the graphical representation includes any one or 
more of the following, Fourier analysis, histogram, signal 
5 moment, surface vibration event analysis and wavelet 
analysis. Other analysis techniques include Sepstrum and 
Homomorphic Deconvolution techniques of non-linear signal 
processing. 

Preferably the method includes the step of 
10 truncating data recorded to an integral number of mill 
rotation periods. 

The method preferably includes the step of 
recording data representing frequency of vibrational events 
occurring over a predetermined period of time. 
15 The method may also include the step of measuring 

or monitoring volumetric load in the machine by identifying 
the toe and shoulder positions of particular matter within 
the machine. 

Preferably the method includes obtaining a polar 
2 0 co-ordinate plot of the volumetric of particulate matter 
within the machine. 

It is preferred that the recorded data is 
processed to provide a display of the location of high 
energy events throughout a polar plot of the machine. 

2 5 It is preferred that the method includes 

recording data for different operational parameters of the 
machine including speed of revolution. 

Preferably the method includes identifying a 
range of angles within which greater 

3 0 numbers /amplitude /frequency of vibrational or other 

mechanical events are occurring. 

It is preferred that a value 

f ' 2k 

35 is able to be produced by the method to provide a 

first approximation for volumetric filling of the machine, 
where _9_ is .the angle (radiance) between the toe and 
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shoulder positions of the particulate material. 

It is preferred that the method includes 
producing a preferred operational speed for the machine in 
order to minimise the number/amplitude/frequency of 
vibrational events occurring within the machine for a 
particular value of 9. 

Preferably the method of controlling operational 
parameters of a machine includes using a value of F for a 
particular amount of particulate matter within a machine in 
order to identify a preferred speed of operation of the 
machine whereby the number/amplitude/frequency of 
vibrational events occurring within the machine is 
minimised. 

Minimisation of number/amplitude/frequency of 
vibrational events occurring within the machine can be 
produced by choosing an operational speed of the machine 
which minimises the value of f . 

It is preferred that impact on the machine caused 
by the toe and shoulder of the particulate matter is able 
to be monitored in order to minimise wear and tear on the 
liner or other parts of the machine. 
Brief Description of the Drawings 

Preferred embodiments of the present invention 
will now be described by way of example with reference to 
the accompanying drawings in which: 

Figure 1 shows a schematic of grinding media 
behaviour in a SAG Mill rotating at abnormally high speed; 

Figure 2 shows a schematic of a basic system for 
monitoring acoustic emissions from a rotating machine; 

Figure 3a shows a graphical representation of an 
accelerometer response on a rotating SAG Mill as a function 
of time; 

Figure 3b shows a plot of standard deviation of 
sampled signal as a function of mill revolution number. 
Manipulated variables are feed rate and mill rotation 
speed; 

Figure 3c shows standard deviation of the sample 
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signal as a function of mill revolution number. 
Manipulated variables are pulp density and ball addition; 

Figure 4a shows a low rotation rate polar contour 
plot of the natural logarithm of a number of vibrational 
• 5 events (E 05 number intervals) as a function of energy 
(volts 2 ) in the radial direction and mill rotation phase 
angle (degrees anti-clockwise from the 3 o'clock position) 
in the azimuthal direction. Contours of large numbers of 
events are at low amplitudes. The mill is rotating 
10 clockwise; 

Figure 4b shows a high rotation rate polar plot 
of the natural logarithm of the number of vibrational 
events (E 0 * 5 number intervals) as a function of energy 
(volts 2 ) in the radial direction and mill rotation phase 
15 angle (degrees ant i -clockwise from the 3 o'clock position) 
in the azimuthal direction. Contours of large numbers of 
events are at low amplitudes. The mill is rotating 
clockwise; 

Figure 4c shows a diagram of amplitude weighted 
20 average event phase angle as a function of revolution 
number and operating conditions; 

Figure 5 shows a graphical representation of 
power spectral density plot upto 1,000 Hz, the parameters 
are 2 16 FFT length, 2 15 number of samples overlap, FFT length 
25 Banning windowing and nil detrending; 

Figure 6 shows a graphical representation of 
power spectral density plot upto Nyquist frequency, the 
parameters are 2 10 FFT length, 2 9 number of samples overlap, 
FFT length Hanning windowing and nil detrending; 
30 Figure 7 a shows a histogram of the natural 

logarithm of the total number of vibrational events 
occurring within an SAG Mill, as a function of event 
amplitude; 

Figure 7b shows a decimated signal for 18 
35 revolutions of data. Decimation of raw signal by a factor 
of 10 3 by successive application of factor of 10 re- 
sampling and filtering; 
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Figure 7c shows a PSD plot upto 1Hz from 
decimated signal (2 11 FFT length and Hanning windowing, 2 10 
number of samples overlap, and nil de trending . ) ; 

Figure 7d shows a spectrogram plot upto 25,000 Hz 
5 (2 12 FFT length and Hanning windowing, 2 11 number of samples 
overlap and nil detrending); 

Figure 7e shows a spectrogram plot upto 1,000 Hz 
(2 12 FFT length and Hanning windowing, 2 11 number of samples 
overlap, and nil detrending) ; 
10 Figure 8 shows a graphical representation of a 

SAG Mill gross power as a function of the standard 
deviation of surfaces for three different mill speed 
settings; and 

Figure 9 shows a graphical representation of pulp 
15 density (percent solids) as a function of the standard 
deviation of surface vibrations for two different sets of 
mill conditions* 

Preferred Embodiment of the Invention 

An example of a system for monitoring vibrational 
20 events occurring within an SAG 10 Mill will now be 
described. 

An SAG Mill 10 consists of a metal cylindrical 

drum containing metallic ball bearings and particulate 

matter to be crushed. 
25 The surface of the mill 10 is provided with an 

accelerometer 11 which is movably fixed thereto and 

includes a radio transmitter. 

A solar panel 12 is electrically connected to the 

accelerometer 11 to provide a power supply. 
30 A radio receiver 13 is connected to a stationery 

part of the SAG Mill framework and receives data 

transmitted by the radio transmitter of the accelerometer 

11. 

The radio receiver is hardwired to a terminal 
3 5 block 14 and a laptop computer 15 is able to be connected 
to the terminal block to receive data sensed by the 
accelerometer. 
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One or more proximity switches are located on the 
rotating part of the mill 1:0 to enable the position of the 
accelerometer to be located relative to the rate of 
rotation of the mill. 

■5 The operation of a SAG mill 10 results in the 

generation of high frequency surface waves on the outside 
of the rotating shell due to collision events within the 
mill. Monitoring of surface waves with a sensitive 
accelerometer therefore provides information on events 

10 inside the mill, particularly impacts of grinding media on 
the liner The grinding media is defined here as non-ore 
(often steel balls) and ore particles (the larger size 
fractions of which contribute to impact grinding) . However, 
measurements of surface vibrations on the outside of the 

15 shell do not simply refledt local impact events on the 
inside of the liner. All of the components of the mill 
behave to some extent as elastic media, permitting the 
propagation of waves generated by collision events within 
the mill, \ flexing' of the mill shell during rotation and 

2 0 external sources such as the drive motor and girth gear. 

Transverse surface vibrational waves propagate around the 
inside of the liner and-around the outside of the shell. 
Meanwhile longitudinal sound waves travel through the 
charge and between the liner and shell. An accelerometer 
25 mounted on. the outside of the shell registers normal 
acceleration due to waves transversely propagating around 
the shell. These waves are damped in accordance with the 
properties of the elastic media between the point of wave 
registration and the origin of the causative event. Hence 

3 0 vibrational events as measured by an accelerometer can be 

expected to be due to causative events over a limited range 
of locations within the mill and associated assembly. 
However, for preliminary analysis it is assumed that the 
vibrations are locally, generated by collision events 
3 5 inside the mill, adjacent to the accelerometer. 

The behaviour of grinding media in a rotating SAG 
mill 10 is usually characterised in a similar manner to 
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balls in a ball mill. Grinding regions are expected to 
consist of the following (see Fig. I): 

• shearing layers of media hear the base of the 
charge article shatter and cleavage (region 1) . 

5 • tumbling media from high to low gravitational 

potential near the top of the charge causing particle 
shatter (region 2). 

• cataracting media impacting in the 'toe' region 
of the. charge causing particle shatter (region 3). 

10 • some abrasion breakage between the grinding 

media and the liner in the contact region between the 
charge and the liner between the w toe and shoulder" (nil 
net force) of the charge (region 4) . 

Higher mill rotational speeds cause cataracting 

15 grinding media to directly impact on the liner above the 
•toe' of the charge (Fig. 1, section 5). Liner 'Wear is 
approximately proportional to the square of the mill speed 
and increases with grinding media size due to the 
consequent increase in impact energy. It is known that the 

20 presence of mineral slurry is very effective in damping the 
impact force of grinding 

media on the liner. Hence liner cracking is interpreted as 
evidence of a problem of excessive direct impact of 
grinding media on the liner. The lifter profile also plays 

25 a strong role in determining the trajectory of cataracting 
grinding media and hence both the location and energy of 
impact of the grinding media on the charge or the liner 
wall. A rectangular lifter profile is thought to result in 
the widest profile of cataracting events for a given mill 

3 0 speed and charge volume This could be interpreted as 
meaning that a rectangular lifter profile leads to the 
highest inherent likelihood of large impact events directly 
on the liner. The composition of the liner itself directly 
influences both impact and abrasion resistance. Use of a 

3 5 liner material that is strongly resistant to abrasion often 
results in low impact resistance and hence an increased 
propensity to cracking. Other causes of liner wear are 



WO 01/03840 



PCT/AU00/00821 



corrosion (in wet mills) and abrasion. However, these sorts 
of wear are not expected to be manifested as liner 
cracking. A non^ intrusive means of quantifying the spatial 
position and intensity of the various types of grinding 
behaviour in SAG mills would be very useful for process 
monitoring and control. SAG mill operators are very keen to 
use a technique that provides a reliable measure of mill 
load Monitoring of the frequency and spatial position of 
large energy particle impacts would also be very useful for 
monitoring and controlling SAG mill liner wear. 
Apparatus 

The surface vibration monitoring system can be 
configured in two ways - basic or advanced. The following 
provides a description of each configuration. 
Basic system 

The surface vibration monitoring apparatus 
consists of a Bruel and Kjaer accelerometer type 4393 
connected to a Bruel said Kjaer charge amplifier type 2634. 
The output from the charge amplifier is connected to a 
microphone beltpack transmitter (AKG type BT5 1) powered by 
a 12 volt gel-cell rechargeable battery. Battery recharging 
is achieved using two solar panels mounted on opposite 
sides of the mill. Transmitted data is received using a 
microphone wireless receiver (AKG type SR5 1) with two 
modified extended antennae. 

Receiver output is connected to a terminal block 
and ribbon cable. The ribbon cable is connected to fast 
data acquisition PCMCIA DAQCard-Al-1 16E -4 linked to a 
laptop computer. 

A magnetic proximity pad was mounted on the mill 
;at 3 o'clock looking from the discharge end. The 
detector/switch was mounted off the mill and connected to 
the terminal block mentioned above. The switching signal 
from the proximity detector is used as a trigger for 
logging of the accelerometer signal. 

The software (written in LabView) can be 
triggered manually or digitally. Triggering occurs when the 
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magnetic pad and detector/ switch comes in close proximity 
as the mill rotates and a 5 volt signal changes to 0 volts. 
Data acquisition then begins with the data being read into 
a rolling buffer and stored into five data arrays. 
Acquisition rate is adjustable up to 100k samples/s. At the 
end of the acquisition process, the mean of the arrays is 
calculated as well as the standard deviation. This; data is 
stored in a file retrieved using Excel and contains the 
following information: 

• Time in seconds from when the system first 
began logging. 

• Signal mean in volts. 

• Signal standard deviation in volts 

• Sample size. 

Unprocessed raw data was also saved in binary 
format for further data processing analysis as mentioned 
later* 

Advanced system 

The advanced system includes the basic system 
described above plus one or more additional accelerometers 
and associated equipment . The accelerometers are of the 
same type as described above. The main difference with the 
advanced system is an improved radio transmission and 
reception apparatus to enable data acquisition down to 
lower frequency ranges (around 1 Hz) . The apparatus 
consists of an Adam module digital radio modem and 
reception system linked to an on-mill computer that 
performs preliminary signal processing before transmission 
to the logging computer. The improved system can run in 
parallel beside the basic system and data is logged and 
saved on the same laptop. The advanced system requires 
extra power for operation so extra solar panels are 
required compared with the basic system. The advanced 
system can consist of multiaccelerometers either in basic 
or advanced configuration to investigate low frequency 
events, event spatial localisation, and events occurring on 
the ends of the mill. Signal data processing techniques are 
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then used to determine surf ace vibration features for given 
operating conditions. A schematic of the basic system 
attached to the shell of a SAG mill is shown in Figure 2. 
EXPERIMENTAL DESIGN 
5 Two series of surface vibrations monitoring runs 

were conducted on the SAG mill at the Red Dome gold mine, 
using the basic system of apparatus. The first was a 
limited series of test runs at different mill operating 
conditions essentially to test the system and establish 

10 that under severe plant duty the system produced data that 
could be processed. Qualitative evaluation of the data 
clearly indicate that features of the data (in the form of 
a voltage-time trace) changed at different operating 
conditions. Features identified were: the minimum signal 

15 strength, length of time between start of signal and 
proximity signal, length of time for events to return to 
background, maximum signal size, number and position of 
high frequency /high energy spikes, amplitude of high energy 
spikes and variation between successive revolutions. These 

20 results highlighted the necessity for further surface 
vibrations monitoring of the Red Dome SAG mill using a 
conditional experimental design so that the influence of 
only one variable could be measured in the context of 
changes in surface vibration features. Thus, this approach 

2 5 would better characterise the potential relationships 

between features and operating settings. 

The second series of test runs at Red Dome 
conformed to the conditional experimental design approach. 
A total of 23 test runs were conducted to investigate how 

3 0 surface vibration features changed with one manipulated 

operating variable at a time. The manipulated operating 
variables were tonnage rate, mill speed, mill discharge 
density and ball addition. In addition to acquiring surface 
vibration information, both control system data and 
3 5 physical plant measurements were taken at each set of 
conditions to confirm test run validity. 

Faster acquisition speeds were possible for the 
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second series of tests and a range of data was collected at 
each manipulated mill condition, A 'normal/ run was 
conducted first at a scan rate of 5 x 10 4 scans/s with raw 
data and statistical averages of the mean and standard 
deviation being saved. Then a 'fast/ run f with acquisition 
rate increased to 1 x 10 5 scans/s, was conducted for a 
duration of 1 X 10 6 scans. Following the fast run, a run 
using the original Labview VI written for the first series 
of lest runs was conducted at the same speed as a 'normal' 
run listed above and saved to a binary file. Lastly, a 
'long' run was conducted with the acquisition speed set at 
5 x 10* scans/s for a duration of 1 X 10 6 scans. Listed 
below are the range of conditions for the manipulated 
operating variables: 

Manipulated variable - tonnage (range 17 0-200 tph) 
Mill speed - 11.8 rpm 
No ball addition 

• Pulp density - 72% solids w/w. 
Manipulated variable - speed (range 12.3-13.8 rpm) 

Tonnage - 2 10 tph 

• No ball addition 

• Pulp density - 72% solids w/w. 
Manipulated variable - ball addition (no-yes 1 kibble) 

• Tonnage - 2 10 tph 

• Mill speed - 13.8 rpm 

• Pulp density - 72% solids w/w. 
Manipulated variable - pulp density (65-72% solids w/w) 

• Tonnage - 210 tph 
Mill speed - 12.5 rpm 
No ball addition 

Manipulated variable - pulp density (72-825$ solids w/w) 

• Tonnage - 210 tph 

• Mill speed - 14.5 rpm 

• No ball addition 

A number of other tests were conducted at 
intermediate conditions for tonnage and speed within the 
range listed above. 
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DATA ANALYSIS TECHNIQUES . 

iThe goal of the data analysis techniques is to 
derive quantitative measures and qualitative visualisations 
based on the response of the accelerometer to shell 
5 vibrations that can be correlated with SAG mill operating 
conditions. This enables vibrsLtion measurements to be used 
for process condition monitoring and as an input to unit 
control. The measures are also useful for inference of the 
rate of liner wear as a function of operating conditions in 

10 with the SAG mill. 

Surface vibrational waves as registered by an 
accelerometer are characterised by a wide variety of 
measures. The first step in data processing prior to 
deriving any of these measures is to truncate the data to 

15 an integral number of mill rotation periods. This is done 
in order to ensure that there is no bias in the data due to 
the sensor detecting changes in mill conditions as a 
function of rotational position of the outer shell. 

The concept of a shell surface vibration event is 

20 important in the data processing. Such an event is defined 
as a positive deviation from nil accelerometer response. 
The amplitude is taken as the peak accelerometer response 
associated with a positive acceleration. This is in 
accordance with a propagating surface wave inducing a 

25 positive acceleration in an accelerometer corresponding to 
a normal stress outwards from the shell. It is hypothesised 
that collision events within the mill, particularly 
grinding media/liner events, will induce a strain that will 
propagate as a wave to the outside of the shell and be 

30 initially sensed as a positive acceleration. A wave train 
due to a collisional event should be composed of an initial 
relatively large, positive acceleration followed by 
negative and positive oscillations of rapidly decreasing 
amplitude. It is expected that the accelerometer will 

35 detect only the first few oscillations of any wave train 
associated with a particular collisional event. Negative 
accelerations are interpreted as part of a wave train 
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belonging to . a previous positive acceleration ^4 are bence 
discarded in terms of registering distinct events. 
Subsequent positive oscillations in a wave train are 
expected to be highly damped due to the low elasticity and 
♦ 5 high damping properties of the liner and the out^r shell. 
Hence it is reasonable that each sequence of positive 
acceleration defines a vibrational event caused by a 
particular media/media or media/ liner collision within the 
mill. 

10 The various measures used to characterise surface 

vibrational waves are as follows: 

1 Mean and standard deviation of the sampled 



signal. 
15 signal, 



2. Power spectral density of the sampled 



3. Histograms of sampled signal amplitude. 

This includes histograms of sampled signal, absolute value 
of sampled signal and the natural logarithm of the absolute 
value of th^ sampled signal. 
20 4. Total number of signal samples and the ratio 

of large to small amplitude samples. The cut-off amplitude 
between large- and small-scale accelerometer responses is 
user defined. 

5. Mean and standard deviation of the amplitude 
25 of surface vibrational events. 

6. Mean and standard deviation of the phase (in 
terms of the position of the accelerometer in the rotation 
cycle of the mill) of surface vibrational events. 

7. Mean and standard deviation of the phase (in 
30 terms of the position of the accelerometer in the rotation 

cycle of the mill) weighted by amplitude of surface 
vibrational events. 

8. Histograms of surface vibrational event 
magnitude. This includes histograms of amplitude and the 

3 5 natural logarithm of the amplitude of events 

9. Total number of events and the ratio of 
large to small amplitude events. The cut-off amplitude 
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between large- and small-scale accelerometer responses is 
again user defined. 

10. The energy of the sampled signal as derived 
from the power spectral density, in the frequency bands 0- 

5 100 Hz, 100-300 Hz and 500-700 Hz and around 18 kHz. These 
frequency bands are deemed by experience to contain most of 
the surface vibrational wave energy information that varies 
with mill operating conditions. 

11. The total energy of the sampled signal 
10 derived from both the power spectral density and the 

amplitude versus time accelerometer response. 

12 . Contour plot of vibrational event numbers as 
a function of SAG mill phase angle and event amplitude. 
Event amplitude and phase angle are identified as 

15 associated with the maximum positive excursion of the 
accelerometer response for any particular event. The 
natural logarithm of surface vibrational event numbers is 
plotted in order to emphasise the position angle of 
relatively large amplitude but infrequent events. 

20 13. Spectrogram time -dependent frequency 

analysis of the sampled signal. 

14. Decimated (decreased sampling rate) versions 
of the signal and associated power spectral density 
analysis. 

2 5 All the above measures are derived for each 

revolution of a continuous monitoring period and over the 
entire integral number of mill revolutions of the same. The 
signal analysis software has been implemented in the MATLAB 
technical computing language. However, any and all of the 

3 0 components of the software could be implemented in -a 

variety of other programming languages. 

A typical accelerometer response trace as a 
function of lime for a rotating SAG mill is shown below 
(Fig. 3) . In this case four full mill revolutions of data 
35 (5 x 10 4 samples per second), previously defined as a long 
run, are recorded. There is clear evidence of periodicity 
in the amplitude of events registered by a single 
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accelerometer as a function of mill rotation angle. 

The standard deviation of the surface vibration 
sampled signal is an important measure of the condition of 
the material within the moving machine. Figures 3b and 3c 

' 5 . show the standard deviation of the sampled signal from a 
SAG mill over a single revolution as a function of mill 
revolution number for a range of manipulated operating 
variables. Figure 3b shows that on average, the signal 
standard deviation is substantially higher at the low 

10 rather than the high feed rate and at the high rather than 
the low mill speed. The former result is thought to 
indicate a low filling level under low feed rate conditions 
at dynamic steady state and hence relatively poor 
cushioning of high energy grinding media impact events on 

15 the liner. The latter is thought to be due to the increase 
in both energy and frequency of grinding media impacting 
directly on the liner. 

Figure 3c shows that the addition of balls 
actually decreases the variability of the standard 

20 deviation across mill rotation periods, but increases its 
mean value. Figure 3c also shows that at high speed, a high 
pulp density condition leads to a decreased standard 
deviation. This may be due to increased dampening of 
grinding media collisions with pulp density. At low speeds 

25 the standard deviation actually increases with pulp 
density. Both Figures 3b and 3c demonstrate considerable 
quasi-periodic variability of the signal standard deviation 
between rotation periods. This may be due to non-steady 
motion of the total charge at frequencies less than or 

30 equal to the rotation of the mill, even when dynamic steady 
state power draft conditions apply. 

A polar contour plot of vibrational event numbers 
as a function of SAG mill rotation phase angle and event 
magnitude allows vibration events to be identified with 

3 5 particular locations in the rotation cycle. Differences may 
be identified in SAG mill operation both between rotation 
periods and with changes in mill operating conditions. The 
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position of either the eyent amplitude or energy weighted 
average phase angle of acoustic ey<ents is a simple and 
valuable quantitative measure in this regard. 

Figure 4a shows: a surface vibration (loosely 

' 5 termed acoustic emission or AE) even|. numbers polar contour 
plot for an "average" SAG mill revolution, based on 19 
single revolutions of data obtained under low speed 
experimental conditions . There is clear evidence of the 
expected localisation of large events in regions where the 

10 charge is thought to be in contact with the liner. In 
particular, there are greater numbers of relatively high- 
energy events in the -290-330 degrees region where it is 
expected that cataracting media impact on the charge. These 
events are thought to identify the position of the "toe" of 

15 the charge. Hence it is inferred that there is strong 
damping of AE waves as they propagate around the shell. 
However, Figure 4a also shows that there is a registration 
of lower strength AE events in regions where the charge and 
grinding media are not expected to be in contact with the 

20 shell liner. These lower energy signals are most likely due 
to surface AE waves propagating around the shell from other 
regions. The prominence of these signals at positions above 
that expected for the n toe" of the charge is in accordance 
with the view that as the sensor rotates through these 

25 positions it is actually registering events that have 
originated further down the shell. Nevertheless, it seems 
likely that very high-energy events recorded by the 
accelerometer at a particular phase angle do reflect 
collisions in the adjacent region of the inner liner. The 

3 0 identification of the boundaries of contact of the SAG mill 
charge with the liner is potentially important as this 
information may be used to deduce the charge volume and 
hence indicate changes in internal conditions as changes in 
operating conditions occur. 

35 Figure 4b is an AE event numbers polar contour 

plot for high rotation rate conditions (14.5 rpm, 12% pulp 
density and- 210 t/h) . . The_- high- energy . events have now 
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clearly split into two peaks at -290 degrees and -330 
degrees (this is even clearer for plots based on individual 
rotation periods). These peaks are postulated to represent 
respectively the *toe' of the charge and a region where 
5 cataracting grinding media directly impinge on the liner. 
The monitoring technique can therefore discriminate large 
energy impacts above the n toe" of the charge. This 
detection of grinding media direct impacts on the liner can 
be used as a predictor of liner wear rate. 

10 Figure 4c shows the amplitude weighted average AE 

event phase angle over a revolution as a function of 
revolution number. As expected from examining Figures 4a 
and 4b, the average event phase angle is in the quadrant 
associated with the "toe" of the charge. There are clear 

15 differences between the average phase angle for the 
different operating conditions previously mentioned. The 
average position angle associated with low feed rate 
conditions is substantially less than the same measure for 
high feed rates. This is physically reasonable, as one 

20 would expect the steady state volumetric loading to be less 
under low feed rate conditions. The average phase angle for 
a high rotation rate is clearly larger than the 
corresponding measure for a low rotation rate. Again, this 
is physically reasonable as one would expect a high 

25 rotation rate to result in more grinding media impacting 
higher up the liner wall. Increasing the pulp density or 
adding grinding balls while maintaining other operating 
conditions lowers the average position angle. These are 
plausible results if increases in pulp density and ball 

3 0 addition are considered to stabilise the bulk behaviour of 
the charge. The significant changes in average event phase 
angle that occur with mill rotation number may indicate 
bulk movement of the charge. 

Figure 5 shows a typical power spectral density 

35 plot obtained for surface vibration waves detected by an 
accelerometer, based on the same set of data. Spectral 
features are readily apparent near the DC channel (< 100 
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Hz) , around 100-200, Hz and near 600 Hz. Finer spikes in 
the power 'spectrum can be seen at frequencies belpw , about 
400 Hz. Figure 4d shows P5p plots for extremes in mill 
speed and feed rate. Spectral features that are sensitive 
5 tp mill operating conditions are apparent , for frequencies < 
100 Hz, around 100-200 Hz and near 600 Hz. Surface 
vibration power is higher at low feed rates at frequencies 
> -100 Hz. High mill speeds result in increased surface 
vibration power at low frequencies (<-100 Hz). The total 

10 power associated with the signal PSD is substantially 
higher for both the low feed rate (-40%) and high speed 
(-30%) operating conditions in comparison with values for 
the respective high f eed ; rate and low speed operating 
conditions . Sharper spikes in the power spectrum can be 

15 seen at frequencies below -500 Hz in the case of high mill 
rotation speed. All these features are well above the 
background noise level and are probably related to bulk 
motion of the charge. In Figure 5, the frequency range was 
restricted vto a maximum of 1000 Hz, because -80% of the AE 

20 signal power was within this range. Power spectral density 
plots have also been obtained up to the Nyqvist frequency 
(half the sampling rate). 

Fig. 6 shows an example based on the same data. 
A prominent spectral feature is apparent at a relatively 

25 high frequency (around 18000 Hz) and experience has shown 
that this feature is sensitive to mill operating 
conditions. Fourier analysis is used in order to try to 
represent the accelerometer response to surface vibrational 
waves in terms of a superposition of sinusoidal waves at 

30 characteristic frequencies. As vibrational events are 
discrete entities, another fruitful approach might be 
wavelet analysis. 

Figure 7a shows a histogram of the number of 
vibrational events as a function of the amplitude of the 

3 5 accelerometer response associated with the event. This is 
an important way of demonstrating both the total number and 
-distribution with amplitude of vibrational events as a 
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function of mill operating conditions. In this manner 
operating conditions that lead to a relatively large number 
of very high amplitude vibrational events can be easily 
identified. Such cases correspond to conditions of high 
5 liner wear. Conversely/ concentration of events at a 
relatively low amplitude indicate ineffective particle 
grinding within the SAG mill. 

Power spectral density analysis associated with 
multiple revolutions of surface vibration data can be used 

10 to obtain information on periodicities in charge behaviour 
at very low frequencies - even down to and below the 
frequency of rotation of the mill. The process used for 
this type of analysis is successive decimation of the raw 
sampled signal in order to decrease the sampling rate of 

15 the data. The decimation process filters the input data 
with a lowpass filter and then resamples the resulting 
smoothed signal at a lower rate. Figures 7b and 7c show the 
decimated signal and power spectral density plot obtained 
via decimation of eighteen revolutions of the surface 

20 vibration signal prior to the addition of grinding media. 
Successive decimation of the signal by factors of ten was 
used to reveal energy peaks at the rotation frequency of 
the mill (0.23 Hz) and several higher harmonics. Spectral 
features sensitive to mill operating conditions are present 

25 over the entire range of frequencies up to those shown in 
Figs. 5 and 6. 

Time-dependent frequency analysis can be used to 
gain further knowledge of the internal behaviour of the 
mill. A spectrogram computes the windowed discrete-time 

3 0 Fourier transform of a signal using sliding windowing. 
Figures 7d and 7e are spectrograms of a single rotation 
period (~5 seconds) of data associated with the low feed 
rate mill operating condition. The colour intensity plot 
displays the power of the signal at each sampling frequency 

3 5 and time. 

Figure 7d shows strong spectral features below ~5 
—kHz- and a distinct feature at -18kHz. This is in accordance 



WO 01/03840 



PCT/AU00/00821 



with power spectral density analysis as demonstrated in 
Figure 6. However, such plots also show that the frequency 
distribution of the power of the signal is not the same at 
ail times over a rotation period. This is more clearly 

• 5 demonstrated in Figure 7e, which shows that power in the 
range -200-1000 Hz is largely restricted to a distinct time 
range in the mill rotation period. The spectral feature at 
-600 Hz previously identified by power spectral density 
analysis can therefore be used as a measure of the 

10 boundaries of the charge during the mill rotation period, 
corresponding to both the 'shoulder' and 'toe' position of 
the charge. Analysis of spectrogram (and power spectral 
density) plots over many rotation periods and for different 
mill operating conditions reveals the sensitivity of these 

15 measures to operating conditions. Spectrogram plots can 
also be represented in a polar fashion (similar to Figures 
4a and 4b) in order to directly determine the shoulder and 
toe positions in terms of phase angles. 

The boundaries of the SAG mill charge as 

20 determination of the weighted average surface vibration 
event phase angle can then be used to provide a first 
approximation for the volumetric filling of the mill. 
Simple geometrical considerations lead to the following 

(fl-sinfl) 
IK 

expression for the fractional filling / of the mill. 
25 Here 6 is the angle (radians) between the toe 

and shoulder positions of the charge. 

CORRELATIONS BETWEEN SURFACE VIBRATION MEASURES AND MILL 
OPERATING CONDITIONS 

As mentioned in the previous section, mean and 

3 0 standard deviation of sampled signal was calculated for all 
test runs. It was postulated that the standard deviation 
would be a useful measure of activity in the mill at a 
given set of conditions. The traces show that the amplitude 
of the signals and the number of high amplitude (energy) 

35 signals changes for different conditions. These changes are 
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then likely to change the standard deviation thereby 
creating an opportunity to establish relationships between 
the surface vibration features encompassed in the standard 
deviation and actual mill operating variables. 
5 8 shows a plot of standard deviation and mill 

gross power for three different mill speed settings. 

Under normal operating conditions the mill gross 
power corresponds to the overall load level in the mill i.e 
the higher the gross power, the higher the load level. 

10 Other manipulated variables, namely ball addition and pulp 
density, were held constant. At each speed, the 
relationship between standard deviation and gross power is 
linear with a negative slope. An increase in the mill gross 
power (load level) leads to a decrease in the standard 

15 deviation of the signal. This result is consistent with 
experience where the mechanism at play is thought to be 
increased damping as a result of the higher load level in 
the mill. The fitted ecjuation for each speed differs and is 
likely to be indicative of some other effects that are at 

20 work. Similar relationships are being developed for other 
surface vibration measures listed in the previous section. 
The relationships shown in Fig. 8 mean that for a given 
speed, the standard deviation measured relates directly to 
the gross power of the mill which itself is an indirect 

2 5 measure of the load level . 

Figure 9 shows a plot of the standard deviation 
against pulp density at two different mill conditions. 

As expected, the relationships for the two 
conditions differ but are consistent with experience. 

3 0 Usually, at higher densities, the thicker slurry inside the 

mill acts as a more effective damper thereby reducing the 
severity of impacts on the shell. This assertion is 
reflected in the measured readings, which show lower 
standard deviation at higher density. Similarly at lower 
3 5 densities, higher surface vibration standard deviation was 
measured . 

ALTERNATIVE SIGNAL ANALYSIS AND CORRELATION DETERMINATION 



WO 01/03840 PCT/AUOO/00821 
TECHNIQUES 

The list of signal analysis techniques used above 
to characterise the surface vibration signals ( including 
techniques 1 based on Fourier analysis, r histogram, signal 
moment and surface vibration event analysis) is by no means 
exhaustive. Wavelet analysis could also be productive given 
that the surface vibration signals are believed to be 
largely due to discrete (time bounded) collision events 
within the mill. The cestruin analysis and homomorphic 
deconvolution techniques of non-linear signal processing 
could also be utilised to further analyse the signals. 
Analysis of signals from multiple sensing devices at 
different locations on a single moving machine could entail 
the use of cross correlation" / cross spectral density 
analysis in order to identify the position of origin (and 
initial intensity) of a vibration. 

Correlations between features of the signal and 
the mill operating conditions can allow the use of surface 
vibration signals as 'soft-sensors' for machine/process 
performance and estimation of unknown plant variables. 
Correlations determined by regression-based analysis of 
surface vibration signals could be extended to include 
consideration of relationships between operating conditions 
and spectral features determined by various types of 
Fourier and/ or wavelet analysis, additional moments and 
statistical measures of the sampled signal and features 
determined by event and histogram statistical analysis. 
Multiple regression and principal component analysis could 
be used to further investigate these linkages. Other 
intelligent analysis methods such as neural networks, 
genetic algorithms, self -organising method, fuzzy logic, 
cluster analysis, Kalman Filter, expert system and 
ARMAX/NARMAX regression-based models could be used to 
analyse the data and seek correlations relating to both 
operating conditions, other process features such as the 
charge particle size distribution and ultimately be used 
for process optimisation and control. 
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Linkages could also be sought to discrete element 
models and/or other models including phenomenal logical 
knowledge. 

ALTERNATIVE TECHNOLOGIES AND APPLICATIONS 

Applications of the surface vibration technique 
extend well beyond SAG mills to other grinding mills e.g 
ball mills / stirred ball mills, jet mills etc. The 
technique is also applicable to other comminution equipment 
such as crushers, impactors, and hammer mills. In fact any 
machine that processes material and a requirement exists 
for a better understanding of the mechanisms occurring 
inside both from a processing and condition monitoring 
viewpoint are potential applications of this technique. 
Machines, the operation of which would benefit with the 
application of [his technique, include but are not limited 
to: 

AG/SAG/ball/rod/vibratory mills 

• Gyratory/ cone/ jaw/rolls crushers 

• Vertical shaft impactors 

• Hammer mills 

• Vertical spindle mills 

• Hydrocycl ones /dense media cyclones 

• Spirals 

• Vibrating/DSM/banana screens 

• Vibrating plate separation devices 

• Flotation cells 

• High pressure grinding rolls/ roller presses 

• Any equipment requiring process or condition 
monitoring 

Rotary kilns and dryers and balling drums 
Alternative technologies that may compete with 
this surface vibration technique are thought to be acoustic 
emission sensors, Le microphones, mounted on the mill and 
configured similarly to this technique. This set-up would 
allow acoustic emission features to be correlated against 
events occurring in the mill during each revolution. The 
hardware required to build this type of system, aside from 
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the sensors, is likely to be very similar to the surface 
vibration system described above. Another possibility is to 
enclose a mill with a mounting frame so that the 
microphones are positioned around the entire circumf erence 
5 of the mill but are mounted off the mill. It is somewhat 
doubtful whether this configuration is practicable. 

The arrangements of experimental apparatus and 
signal analysis techniques as applied to .monitoring a SAG 
mill as described in this document have been advanced 

10 merely by way of explanation. Many modifications may be 
made thereto both for further monitoring of SAG mills and 
other material processing machines without departing from 
the spirit and scope of the invention which includes every 
novel feature and combination of novel features herein 

15 disclosed. 
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CLAIMS : 

lb A system for monitoring mechanical waves 
from a machine which in operation has moving particulate 
matter therein, the system including at least one sensor 
located ozv the machine at a location away from the central 
axis of the machine, the sensors being for sensing acoustic 
waves and including a transmitter for transmitting signals 
representing the sensed mechanical waves to a receiver at a 
location remote from the sensor (s), a data processor 
connected to the receiver for receiving signals from the 
receiver which signals represent the mechanical waves and 
processing the signals to produce output signals for 
display on a display means, wherein the output signals for 
display represent one or more parameters indicative of 
mechanical waves emitted from the machine over a 
predetermined period of time. 

2. The system as claimed in claim 1 wherein the 
receiver is located on a stationery surface separate from 
the machine. 

3. The system as claimed in claim 2 including a 
power supply located on the machine. 

4. The system as claimed in claim 3 wherein the 
or each sensor is located on an exterior surface of the 
machine . 

5. The system as claimed in claim 4 wherein the 
data processor is adapted to produce output signals which 
represent a plurality of acoustic events occurring within 
the machine, amplitudes of the acoustic events and data 
relating to the position of the acoustic events. 

6. The system as claimed in claim 4 or claim 5, 
including a plurality of sensors spaced around the 
periphery of the machine to enable polar co-ordinates of 
the origin of emissions to be located. 

7 . The system as claimed in claim 6 wherein the 
sensors are ecpiispaced around the periphery of the machine. 

8. The system as claimed in claim 3 or 4 
wherein sensors are arranged in an array around the machine 
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and along the length of the machine to enable to a three 
dimensional co-ordinate axis to be plotted of the location 
of the origin of omissions from the machine. 

9. The system as claimed in claim 6 or claim 7 
5 including at least one proximity detector for monitoring 

the location of the sensors at a predetermined time, 
whereby data from the proximity switch is adapted to be 
communicated to the data processor. 

10. The system as claimed in claim 7 wherein the 
10 data processor includes a timing means for calculating the 

location of the or each sensor at a predetermined time. 

11. The system as claimed in claim 9 wherein the 
or each sensor includes an accelerometer which is adapted 
to transmit data relating to the frequency of vibrational 

15 events occurring within the machine and the amplitude of 
the vibrational events at particular locations within the 
machine to the transmitter. 

12 . A method of analysing operational parameters 
of a machine having a moving particulate material therein, 

20 the method including the steps of recording data 
representing a number of mechanical events occurring within 
the machine over a predetermined period of time, the 
amplitude of the mechanical events occurring over the 
predetermined period of time and positional data relating 

25 to the position of the mechanical events occurring within 
the machine, displaying a graphical representation of the 
recorded data, the graphical representation including 
parameters relating to the number of mechanical events, the 
amplitude of mechanical events and the position of 

3 0 mechanical events occurring within the machine during the 
machines operation. 

13 . The method as claimed in claim 2 wherein the 
graphical representation of recorded data includes mean and 
standard deviation of vibrational events occurring within 

3 5 the machine, power spectral density of vibrational events 
occurring within the machine and histograms of amplitude of 
vibrational events occurring within the machine. 
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14. The method as claimed in claim 13 including 
the step of measuring volumetric load of the particulate 
matter within the machine by identifying the toe and 
shoulder portions of the particulate matter. 
• 5 15. The method as claimed in claim 14 wherein 

volumetric load is determined from a polar co-ordinate plot 
of events occurring within the machine. 

16. The method as claimed in claim 15 wherein 
volumetric load is calculated for a range of angles in 

10 which events within the machine have greatest deleterious 
effect on the interior of the machine. 

17. The method as claimed in claim 16 wherein a 
value for volumetric filling of the mill is produced from 
the recorded data and the value of volumetric filling 

15 (fl-sinfl) 

f ~ In 

where 0 is the angle (radiance) between the toe and 
shoulder positions of the particulate matter. 

1$. A method of controlling operational 
20 parameters of a machine having a moveable substance 
therein, the method including the steps of recording data 
representing a number of vibrational events occurring 
within a machine over a predetermined period of time, 
amplitude of the vibrational events occurring over the 

2 5 predetermined period of time and position data relating to 

the position of the vibrational events over the 
predetermined period of time, determining zones within the 
machine which are subject to predetermined levels of wear 
and altering the machine operational characteristics to 

3 0 reduce the levels of wear for the zones. 

19. A method of identifying the volumetric load 
of particulate matter within a machine including the steps 
of receiving data, representing a number of mechanical 
events occurring within the machine over a predetermined 
35 period of time, the amplitude of the mechanical events 
occurring over the predetermined period of time and 
positional data relating to the position of the mechanical 
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events occurring within the machine, processing the 
received data to identify toe and shoulder positions of the 
particulate matter within the machine whereby the location 
of maximum deterioration of an inside surface of the 
• 5 machine can be minimised. 

20. The method as claimed in claim 19 wherein 
data is received for a plurality of speeds of the machine. 

21. The method as claimed in claim 20 including 
the step of identifying the fractional filling f of the 

10 machine where 

(g-sinfl) 
; 2/r 

with 8 being the angle (radiance) between the toe and 
shoulder positions of the charge.. 
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